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ABSTRACT 

Strain-mediated muitiferroic composites, i.e., piezoelectric-magnetostrictive heterostructures, hold profound promise for 
energy-efficient computing in beyond Moore’s law era. While reading a bit of information stored in the magnetostrictive 
nanomagnets using a magnetic tunnel junction (MTJ), a material selection issue crops up since magnetostrictive materials in 
general cannot be utilized as the free layer of the MTJ. This is an important issue since we need to achieve a high magnetore¬ 
sistance for technological applications. We show here that magnetically coupling the magnetostrictive nanomagnet and the 
free layer e.g., utilizing the magnetic dipole coupling between them can circumvent this issue. By solving stochastic Landau- 
Lifshitz-Gilbert equation of magnetization dynamics in the presence of room-temperature thermal fluctuations, we show that 
such design can eventually lead to a superior energy-delay product. 


Introduction 

Electric field-induced magnetization switching in strain-coupled muitiferroic composites is a promising mechanism that can 
possibly harness an energy-efficient binary switch replacing the charge-based traditional transistors for our future information 
processing paradigm.'’^ A voltage applied across such devices strains the piezoelectric layer and the generated stress on the 
magnetostrictive layer induces a magnetic anisotropy in it^“'^ and can switch its magnetization.*’^’These straintronic 
devices operate at room-temperature and the study estimates very promising performance metrics, e.g., energy dissipation 
of ^1 attojoule (aJ) and sub-nanosecond switching delay, suitable for technological application purposes.*^’*® Experimental 
efforts to investigate such device functionality has demonstrated the induced stress anisotropy in magnetostrictive nanomag¬ 
nets, *^“^^ while the direct experimental demonstration of switching speed (rather than ferromagnetic resonance experiments to 
get the time-scale) and using low-thickness piezoelectric layers while avoiding considerable degradation of the piezoelectric 
constants [e.g., < 100 nm of lead magnesium niobate-lead titanate (PMN-PT)]^'*’^^ are still under investigation. 

There are proposals on devising both memory^’ *^’ *^ and logic devices using strain-mediated multifeiToic composites 
by energy-efficient writing of a bit of information in the magnetostrictive nanomagnets.However, while electrically reading 
the magnetization state of the magnetostrictive nanomagnet using a magnetic tunnel junction (MTJ),^^“^^ we need to tackle a 
material selection issue since the magnetostrictive materials in general cannot constitute the free layer of an MTJ. The widely- 
used material that is used for the free layer of an MTJ is CoEeB,^® which leads to high tunneling magnetoresistance (TMR) of 
300%.^^ The incorporation of half-metals as the free layer can lead to even better TMR of more than 1000%.^^ To tackle this 
material selection issue, we propose to magnetically couple the magnetostrictive nanomagnet and the free layer of an MTJ, 
e.g., to utilize the magnetic dipole coupling in between them separated by an insulator (see Eig. 1). During write operation, as 
the magnetization of the magnetostrictive layer rotates upon application of stress, the free layer’s magnetization also rotates 
concomitantly and it can be read by an MTJ. Similar methodology of incorporating an insulator for utilizing magnetic dipole 
coupling has been proposed in the context of input-output isolation for logic design purposes.Note that the input-output 
isolation is inherent in muitiferroic devices due to the presence of the insulating piezoelectric layer.^^ We study the effect 
of this dipole coupling by solving stochastic Landau-Lifshitz-Gilbert equation of magnetization dynamics in the presence of 
room-temperature thermal fluctuations. The results reveal that such dipole-coupled design can lead to lowering the energy 
dissipation and a superior energy-delay product. 

Model 

The nanomagnets are modeled in the shape of elliptical cylinders with the cross-sections lying on the y-z plane; the major axis 
points along the z-direction, and the minor axis along the y-direction (see Eig. la). In standard spherical coordinate system, 9 
is the polar angle and (j) is the azimuthal angle. Any deflection out of magnet’s plane (0 = ±90°) is termed as out-of-plane 
excursion. We solve the magnetization dynamics using stochastic Landau-Lifshitz-Gilbert (LEG) equation in the presence of 
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Figure 1. Separating read and write units in multiferroic composite devices, (a) A write unit utilizing multiferroic 
composites, i.e., piezoelectric-magnetostrictive heterostructures, and axis assignment using standard spherical coordinate 
system, (b) A read unit in the form of a magnetic tunnel junction (MTJ) is incorporated to read the magnetization state of the 
magnetostrictive nanomagnet, which acts as the free layer in the MTJ. However, materials like half-metals, CoFeB are 
suitable to constitute the free layer since that leads to a high tunneling magnetoresistance (TMR). Hence, it necessitates to 
decouple the read unit and write unit in a multiferroic device, (c) The read unit and the write unit are separated by an 
insulator but the magnetostrictive nanomagnet and the free layer are magnetically coupled through dipole coupling. Note that 
the read current flows through the read unit (MTJ) only and the detailed contacts are not shown in this schematic diagram, (d) 
A detailed design to make contacts for the read and write units in the proposed design. The read current must flow through 
the read unit (MTJ) only so that the TMR of the MTJ does not get affected. A metallic layer on the insulator layer can be 
incorporated to form an equipotential surface, which makes a contact for the free layer in the MTJ. The write terminals are 
also shown here. Since the piezoelectric layer is much more resistive than that of the insulator, having the insulator in the 
path does not affect the operation and hence it does not necessitate another metallic layer below the insulator layer for the 
corresponding write unit terminal. 
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room-temperature (300 K) thermal fluctuations. Note that there are two nanomagnets (magnetostrictive and free layer) having 
a dipole coupling between them. Stress is generated only on the magnetostrictive nanomagnet, hence we have additional stress 
anisotropy to consider for the magnetostrictive nanomagnet. We will use the subscripts m and / to denote any parameter for 
the magnetostrictive nanomagnet and the free layer nanomagnet, respectively. 

The magnetization Mm (Mf) of the magnetostrictive (free layer) nanomagnet has a constant magnitude but a variable 
direction, so that we can represent it by a vector of unit norm nim = Mm/|Mm| = Cr (nif = Mf/|Mf| = Cr) where Cr is the unit 
vector in the radial direction in spherical coordinate system represented by (r,0,0). The other two unit vectors corresponding 
to the polar angle 9 and the azimuthal angle (j) are eg and e,;,, respectively. 

The potential energies of the magnetostrictive nanomagnet and the free layer nanomagnet can be expressed, respectively, 
as 
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[The dipole coupling between two magnetic moments Mm and Mf separated by a distance vector R can be expressed as^ 

£'d;po;e = (l/4;r^o^^)[(Mm.Mf)-(3//;2)(Mm.R)(Mf.R)], where putting |Mm| |Mf| andR = f;ex, 

we get the equation (5).] 

Ms,m iMs,f) is the saturation magnetization of the magnetostrictive (free layer) nanomagnet, Q.m i^f) is the volume of the 
magnetostrictive (free layer) nanomagnet, Nd-pp,m {Nd-pp,f) is the component of the demagnetization factor for the mag¬ 
netostrictive (free layer) nanomagnet along p-direction, which depends on the nanomagnet’s dimensions,^’(3/2)2^ is the 
magnetostrictive coefficient of the single-domain magnetostrictive nanomagnet,^ a is the stress on the magnetostrictive nano¬ 
magnet, and R is the center-to-center distance between the nanomagnets. 

The initial orientation of the magnetizations is antiparallel due to dipole coupling between the nanomagnets. From equa¬ 
tion (3), note that when we apply a sufficient stress (compressive stress for materials with positive Aj or vice-versa so that 
the product A* O’ is negative) on the magnetostrictive nanomagnet, the induced stress anisotropy can beat the shape anisotropy 
of the nanomagnet and rotate its magnetization toward the hard axis.^’^^’*^ As the magnetization of the magnetostrictive 
nanomagnet rotates due to the stress anisotropy induced in it, the magnetization of the free layer nanomagnet does also rotate 
due to the magnetic dipole coupling between the nanomagnets. The magnetizations keep antiparallel orientation when they 
reach the hard axis. Upon removal of stress from the magnetostrictive nanomagnet, it switches to the opposite direction due to 
out-of-plane excursion of magnetization.The dipole coupling switches the magnetization of the free layer concomitantly. 

The effective held and torque acting on the magnetostrictive nanomagnet due to the gradient of potential landscape as 
given by the equation (1) can be expressed as Hefrm = -'^Etotai,m = -{dE,o,ai,m/d9m)te - {l/sin9,n) {dE,o,ai,m/di)m) and 
Te.ih = Him X Heff.nn respectively. Similarly, the effective held and torque acting on the free layer nanomagnet due to the 
gradient of potential landscape as given by the equation ( 2 ) can be expressed as Hetr,f = —'^Efotai.f = —{dEtotai.f / d9f)eg — 
{\/sin9f) {dE,otai,f/d(pf)ei^ and Ts.r = nif x respectively. 

The thermal held and the corresponding torque acting on the magnetostrictive (free layer) nanomagnet can be written^^’’^*' 
as HTH.m(f) = ee and TxH.mif) = m^{f) x HTH,m(f). respectively, where 

[^U,m(/) yt)s9m{^f) ttOStpm^f^ 4“ ^y,m{f) t^tds9m{^f'^sin(pm{^f^ 5 (6a) 
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is the phenomenological damping parameter of the magnetostrictive (free layer) material, 7 is the gyromagnetic ratio 
for electrons, is the simulation time-step, G„,(^)(0,1) is a Gaussian distribution with zero mean 

and unit variance for the magnetostrictive (free layer) nanomagnet,^* k is the Boltzmann constant, and T is temperature. 

The magnetization dynamics of the magnetostrictive (free layer) nanomagnet under the action of the two torques j„(f) 
and Txjj is described by the stochastic Landau-Lifshitz-Gilbert (LLG) equation^*’’as 
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The magnetization dynamics of the two nanomagnets represented by the equations ( 8 ) and (9) are coupled through the 
dipole coupling [see equation (5)]. These coupled equations are solved numerically to track the trajectories of the two magne¬ 
tizations over time. 

The internal energy dissipation in the magnetostrictive (free layer) nanomagnet due to Gilbert damping can be expressed 
as £d,,„(/) = lo Pd,m{f){t)Bp where t is the switching delay and the instantaneous power dissipation can be calculated as 


Pd,m{f){t) 


«m(/) |y| 


TE.m(f) {t ) 


( 11 ) 


We sum up these two internal energy dissipations Ed,m and Ed,f alongwith the energy dissipation due to applying voltage 
(which is miniscule jq determine the total energy dissipation. 


Results 

The magnetostrictive nanomagnet is made of polycrystalline Terfenol-D and it has the following material properties - sat¬ 
uration magnetization (Ms,m)' 8x10^ A/m, Gilbert damping parameter (am)' 0.1, Young’s modulus (Y): 80 GPa, magne¬ 
tostrictive coefficient ((3/2)Xs): h-90x 10^^ (Refs. 14,44-46), and Poisson’s ratio (v): 0.3 (Ref. 47). The free layer nano¬ 
magnet is made of widely-used CoFeB, which has the following material properties - Gilbert damping parameter («/); 
0.01, saturation magnetization (Ms,/): 8x10^ A/m.^^ The dimensions of both the single-domain nanomagnets are chosen 
as lOOnm x 90 nm x and the center-to-center distance between the nanomagnets is R = 40 nm. For the piezoelectric 
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Figure 2. Magnetization switching in the magnetostrictive and free layer nanomagnets upon application of stress on 
the magnetostrictive nanomagnet, (a) Magnetizations of the magnetostrictive and free layer nanomagnets are 
antifeiTomagnetically coupled due to dipole coupling, (b) Potential landscapes (bold lines correspond to magnets’ planes) 
and positions of the magnetizations for both the nanomagnets during switching. When no stress is applied, the potential 
landscapes of the nanomagnets are monostable due to dipole coupling and the magnetizations are antiferromagnetically 
coupled. If a sufficient stress is exerted on the magnetostrictive nanomagnet, the magnetizations come to their hard-axes and 
still remain antiferromagnetically coupled, however, they are deflected out-of-plane. This out-of-plane excursion eventually 
leads to the full 180° switching when the stress is released/reversed, (c) The LLG simulation results show that both the 
magnetizations successfully switch by 180°. A stress of 30 MPa is applied on the magnetostrictive nanomagnet and the stress 
is reversed to aid the switching speed when the magnetization comes to the hard axis. The ramp (rise and fall) time of stress 
is 60 ps. No thermal fluctuations is considered but the deflection in initial orientations of the magnetizations is taken as 
^3.25°, which is the thermally mean value due to room-temperature (300 K) thermal fluctuations when no stress is active. 
When 9,„ becomes < 5°, the switching is deemed to have completed and the switching delay is recorded as 0.351 ns. 


layer, we use PMN-PT, which has a dielectric constant of 1000, c/3i=-3000 pm/V, and £/32=1000 pm/V.^* We assume the 
piezoelectric layer’s thickness fp,(.20=24 nm (Ref. 14) and thus V = 1.9 mVs (2.9 mVs) of voltages would generate 20 MPa 
(30 MPa) compressive stress [(7 = Yd^fflV/tpiezo), where deff = (c /31 —c/ 32 )/(l -f v)] in the magnetostrictive Terfenol-D 
layer. Note that avoiding considerable degradation of the piezoelectric constants at such low-thickness (24 nm) piezoelectric 
layers is under research.Modeling the piezoelectric layer as a parallel plate capacitor, the capacitance C=2.6 fF and thus 
CV^ energy dissipation turns out to be <0.1 aJ. This is the basis of ultra-low-energy computing using these multiferroic 

devices. *’^44,26,27 

When the magnetizations of the magnetostrictive and free layer nanomagnets are exactly aligned to their easy axes (e.g., 
9m = 180° and 9f = 0°), the torques acting on the magnetizations are exactly zero and hence only thermal fluctuations can 
deflect the magnetizations from their initial orientations. When no stress is active on the magnetostrictive layer, we solve 
the stochastic LLG equation in the presence of room-temperature thermal fluctuations to determine the distributions of the 
magnetizations’ initial orientations and calculate the mean orientations of the magnetizations (^3.25°). The initial distribution 
of magnetization is a Boltzmann distribution and matching the numerically calculated mean orientation of magnetization with 
the one calculated from the equipartition theorem depicts the validity of incorporating thermal fluctuations.^^ 

Figure 2 shows that the magnetizations of the magnetostrictive nanomagnet and the free layer nanomagnet rotate concomi¬ 
tantly upon application of 30 MPa stress on the magnetostrictive nanomagnet. The magnetizations of two nanomagnets come 
to their respective hard axes and remain antiferromagnetically coupled. The initial values of the azimuthal angles and 

are chosen as 270° and 90°, respectively, however, they can be just opposite too, which is equally possible. During 
the course of magnetization dynamics, the exerted stress rotates the magnetization of the magnetostrictive nanomagnet out-of- 
plane and subsequently the magnetization of the free layer also gets deflected out-of-plane due to dipole coupling, as depicted 
in the Fig. 2b. This out-of-plane excursion increases the switching speed tremendously and creates an intrinsic asymmetry to 
facilitate a complete 180° switching of the magnetizations deterministically even in the presence of thermal fluctuations.^^ 
The LLG simulation results as shown in the Fig. 2c depict that both the magnetizations have completed full 180° switching. 

The magnetization switching procedure described above requires to read the magnetization state using MTJ to sense 
when magnetization reaches around 9m = 90° (since room-temperature thermal fluctuations make the traversal time a wide 
distribution), so that stress can be brought down thereafter.*^ Note that there is tolerance around 0,„ = 90°, i.e., stress does not 
need to be withdrawn exactly at 0„, = 90° since it is shown that the internal magnetization dynamics provides such tolerance. *^ 
These are purely dynamical phenomena contrary to steady-state analysis. Any additional element for comparison can be built 
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Case 

Free layer? 

Stress (MPa) 

'^mean 

Xstd (ns) 

Efn (^) 

Ef (aJ) 

£rf(aJ) 

(aJ) 

^total (^J) 

^total '^mean (aj-ns) 

(a) 

No 

20 

0.444 

0.080 

0.89 

NA 

0.89 

0.0297 

0.9296 

0.4127 

(b) 

Yes 

20 

0.529 

0.120 

0.62 

0.07 

0.69 

0.0297 

0.7296 

0.3860 

(c) 

Yes 

30 

0.379 

0.080 

0.81 

0.08 

0.89 

0.0669 

0.9792 

0.3711 

(d) 

No 

30 

0.368 

0.064 

1.09 

NA 

1.09 

0.0669 

1.1792 

0.4340 


Table 1 . Performance metrics for four different cases considered. Cases (a) and (d) [corresponding to 20 MPa and 30 MPa 
stress, respectively] do not consider the additional free layer, while the cases (b) and (c) [corresponding to 20 MPa and 30 
MPa stress, respectively] consider the free layer magnetically coupled to the magnetostrictive nanomagnet. With the 
introduction of the free layer at the same stress level, the switching delay metrics (mean and standard deviation) get worse 
while it dissipates less energy and leads to less Efotai %iiean product. Case (c) has the lowest Efotai %nean- 


using these energy-efficient multiferroic devices in general.^® Note that researchers are trying to replace the traditional switch 
based on charge-based transistors by a new possible “ultra-low-energy” switch (e.g., using multiferroic composites). Therefore, 
any circuitry can be built with the energy-efficient switch itself rather than the conventional transistors. Usually, it requires 
several peripheral circuitry in conjunction with the basic switch in a system.^*’^^ While researchers report on the performance 
metrics of the basic switch itself, the total energy dissipation considering the other required circuitry does not change the order 
of energy dissipation, utilizing the respective devices.This was the understanding while claiming energy-efficiency using 
such magnetization switching mechanism*’^’ and computing methodologiesbased on such switching methodology. 

It may be possible to harness more asymmetry in the system apart from the intrinsic asymmetry due to out-of-plane 
excursion as described above so that the sensing mechanism for dynamic withdrawal of input voltage may not be necessary. 
Interface and exchange coupling can also provide asymmetry during switching,^^ particularly it helps to maintain the direction 
of switching rather than toggling*^ the magnetization direction, and it does not require any sensing procedure. 

To understand the effect of incorporating the dipole-coupled free layer on the performance metrics, we solve the stochastic 
LLG equation^** at room-temperature (300 K) and tabulate the performance metrics in Table 1 for four different cases. For 
each case, we perform a moderately large number (10,000) of simulations and when the magnetization of the magnetostrictive 
nanomagnet reaches 9m < 5°, the switching is deemed to have completed and the switching delay T for that trajectory is 
recorded. Then we determine the following performance metrics: mean value of switching delay (Xmean), standard deviation 
of switching delay {Tstd)^ the mean values of the energy dissipations Em (in the magnetostrictive layer) and Ef (in the free 
layer) due to Gilbert damping in the magnets, the energy dissipation due to applying voltage ‘CV^’ = 3CV^ (since stress is 
reversed'"^’ '^), and the total energy dissipation Efotai =Ed + ‘CV^\ where E^ = Em+Ef. 

For the case (a), we do not have any additional free layer and just consider the switching in the magnetostrictive nanomag¬ 
net, while for the case (b), we do have the free layer. For both the cases (a) and (b), the stress is 20 MPa. We also consider the 
distribution of initial orientation for the case (a) and the mean value of deflection of initial orientation from the easy axis turns 
out to be ^^3.9°. This value is higher than that of the case (b) [^3.25°], when the dipole-coupled free-layer is introduced. The 
reason behind is that the dipole coupling energy confines both the magnetizations more in their respective potential wells so 
the magnetizations’ deflection is less. Hence, while considering switching with dipole coupling, magnetizations on average 
start nearer from the easy axis and therefore it takes more time for switching to be completed. This is reflected in the mean 
value of switching delay if we compare it for the cases (a) and (b). By fixing the same initial orientation of magnetization for 
the cases (a) and (b), it is noticed that dipole coupling in fact speeds up the switching process so the increase in switching 
delay for the case (b) is entirely due to the less deflection in the initial orientation of magnetization as described above. The 
case (b) incurs less energy dissipation in total, which indicates the delay-energy trade-off, i.e., slower switching dissipates less 
energy. Note that case (b) has lower energy-delay product {E,otai '^mean) compared to the case (a). 

To investigate the performance metrics with the incorporation of the dipole-coupled free layer further, we increase the 
stress to 30 MPa and tabulate the results as case (c) in the Table 1 . We note that the mean switching delay has got reduced 
compared to the case (a) while they incur the same amount of energy dissipation due to Gilbert damping. We plot the 
corresponding switching delay distribution for the case (c) in the Fig. 3. Such distribution can be achieved experimentally 
by time-resolved measurements.^'^ It needs to be pointed out that we can generate a maximum amount of stress on the 
magnetostrictive layer dictated by the maximum strain induced in it, so we also consider 30 MPa stress without the free layer 
and tabulate the results as case (d). The mean switching delay Xmean for case (d) is very close to that of case (c), but it has 
the highest energy dissipation Etotai (™d also Etotai'^mean) among the four cases considered, while case (c) has the lowest 
Etotai%nean- Assuming a performance metric Efotai where T = Xniean + 10 Tjrj, the case (c) still would have the lower product 
compared to the case (d). 

Note that the switching delay is optimized for lowest value with respect to dipole coupling strength, which can be tuned by 
varying the thickness of the insulator separating the magnetostrictive nanomagnet and the free layer. With higher thicknesses of 
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Figure 3. An illustrative distribution of switching delay for the case of magnetostrictive nanomagnet. A moderately 
large number (10,000) of simulations performed in the presence of room-temperature (300 K) thermal fluctuations. The 
dipole-coupled free layer is considered in this case and 30 MPa stress is applied with 60 ps ramp (rise and fall) time. This 
wide distribution is caused by the following two reasons: (1) thermal fluctuations make the initial orientation of 
magnetization a distribution, and (2) thermal kicks during the switching makes the switching delay a distribution too; the 
former one has a higher effect than that of the latter. The mean and standard deviation of this distribution are 0.379 ns and 
0.080 ns, respectively. 


the insulator, the magnetizations of the two layers do not quite rotate concomitantly and thus the switching delay is increased, 
while for lower thicknesses, higher dipole coupling rotates the magnetization out-of-plane so much that it leads to precessional 
motion and it increases the net switching delay. 

Discussion 

We have addressed the material selection issue while reading out the state of the magnetostrictive nanomagnet in a multiferroic 
composite. The proposed design provides us the flexibility to use the best materials for the magnetostrictive nanomagnet and 
the free layer in an MTJ separately. Rather than dipole coupling, we can also utilize exchange coupling between the two 
nanomagnets to magnetically couple them. Note that this is a general strategy, which can be also utilized in spin-transfer- 
torque switching of nanomagnets where we can use the switching nanomagnet to be made of CoFeB and the free layer to 
be made of half-metals for higher TMR. Hence, it will motivate experiments and further theoretical studies on this front. 
Moreover, it turns out that this design also enhances the energy-delay performance metric. Such ultra-low-energy and non¬ 
volatile (leading to instant turn-on computer) computing paradigm is particularly promising to become the staple of our future 
information processing systems. 
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